Spin exciton in a quantum dot with spin-orbit coupling at high magnetic field 
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Coulomb interactions of few (N) electrons confined in a disk shaped quantum dot, with a large 
magnetic field B = B* applied in the z— direction (orthogonal to the dot), produce a fully spin po- 
larized ground state. We numerically study the splitting of the levels corresponding to the multiplet 
of total spin S = N/2 (each labeled by a different total angular momentum J z ) in the presence of an 
electric field parallel to B, coupled to S by a Rashba term. We find that the first excited state is a 
spin exciton with a reversed spin at the origin. This is reminiscent of the Quantum Hall Ferromagnet 
at filling one which has the skyrmion-like state as its first excited state. The spin exciton level can 
be tuned with the electric field and infrared radiation can provide energy and angular momentum 
to excite it. 

PACS numbers: PACS numbers:73.21.La,73.23.-b,78.67.Hc 



I. INTRODUCTION 

Quantum dots arc semiconductor devices in which elec- 
trons are confined to a small area within a two dimen- 
sional electron gas by properly biasing metal gates added 
to the structure^. In an isolated quantum dot (QD) the 
confining potential gives rise to quantized single particle 
energy levels. However, electron-electron interaction de- 
termines the dot properties. In dots having a diameter of 
~ lOOnm, the level spacing and the Coulomb energy are 
of the order of ~ lmeV and most charging properties can 
be included within the Hartree-Fock (H — F) approxima- 
tion, just like in atoms. Correlation effects do not signif- 
icantly alter the charging properties, but may strongly 
influence the spin properties of the confined electrons. 
One striking evidence of this is the fact that Hund's rule, 
which is typical of atoms, is often satisfied in dots^. How- 
ever, the reduction of the energy scale by a factor of 10~ 3 
with respect to atoms, enhances the sensitivity of the 
electrons in the dot to an external magnetic field. 
In this paper we discuss the spin properties of an iso- 
lated vertical QD in presence of a magnetic field B in 
the z— direction, orthogonal to the dot disk (cylindrical 
symmetry is assumed). We also include the spin orbit 
coupling (SO) induced by an electric field along z. 

In the cylindrical geometry, orbital effects are 
dominant 4 . Indeed Zeeman spin splitting does not drive 
any spin polarization in these systems and can be of- 
ten ignored 5 . However correlations combine orbital and 
spin effects together and can be probed by magnetocon- 
ductance measurements in a pillar configuration^. Spin 
properties are quite relevant to conductance, in view of 
the possibility of spin blockaded, Kondo effec1pi2iii, or 
Berry phase induced tuning 12 . 

Quantum numbers labeling the dot energy levels are 
the number of electrons N, the total orbital angu- 
lar momentum along z, M, the total spin S and the 



z— component of the spin S z . By increasing the magnetic 
field B, both M and S increase. Finally, at B = B*, a 
fully spin polarized (FSP) state is reached. The increas- 
ing of the total spin S was measured in a dot with about 
30 electrons, a striking evidence of e-e correlations^. 
While, in the absence of interactions, the density of the 
FSP state becomes uniform over the dot area, which is 
contracted to a minimum, in the real case, the e-e inter- 
action tends to reduce the density at the center of the 
disk, by compressing the electrons at the dot edge (see 
FigEl [middle panel]). By further increasing B, the elec- 
tron density reaches a maximum value. The dot becomes 
a so called "Maximum Density Droplet" (MDD)2±. For 
larger B values, the FSP state is disrupted: the dot den- 
sity reconstructs i.e. an annular local maximum of the 
density is produced at the edge of the dolii 5 - with break- 
ing of azimuthal symmetry at the edge (de Chamon-Wen 
phased) . 

Various numerical calculations 1 ^ have investigated 
these subsequent electronic transformations which ap- 
pear as crossing of levels with different quantum num- 
bers. H — F calculations are known to incorrectly favor 
spin-polarized states**. Spin density functional calcu- 
lations have been performed for dots including a larger 
number of electrons 1 ^. The density functional approach, 
with a good choice of the parameters of the poten- 
tials, can reach a significative agreement with the exper- 
iments, but it may introduce uncontrolled approxima- 
tions. When the electron density is reduced, a Wigner 
molecule can be formed. Recently, this broken symme- 
try state has been studied in the absence of an external 
magnetic field, using Quantum Monte Carlo simulations, 
with a multilevel blocking algorithm which is free of the 
sign problemSS. 

In this work, we use exact diagonalization for few elec- 
trons with azimuthal symmetry"? to discuss the spectral 
properties of the FSP dot which is stabilized by the 
Coulomb interaction. Electrons are confined to a two- 
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dimensional (2— d) disk by a 2 — d parabolic potential and 
interact via the full Coulomb repulsion whose strength is 
parametrized byU = e 2 / kI. Here I is the magnetic length 
due to the parabolic confinement in presence of a field B 
along z and k is the static dielectric constant. The con- 
finement of the electrons in the x — y plane implies the 
presence of an electric field in the z— direction, provided 
by the band bending of the hetherostructure. This gives 
rise to the so called SO Rashba termSi, which can be 
enhanced even further in a non-linear conductance mea- 
surement, when an extra bias voltage V s d is applied to 
the contacts of a vertical structure. 

In the presence of SO, J z = M + S z becomes the 
good quantum number. The FSP ground state (GS) 
has J z = N(N -l)/2 + S (J z = 25/2 for N = 5), while 
the first excited state (denoted as SKD state in the fol- 
lowing) has J z = N(N - l)/2 + S - 1 (J z = 23/2 for 
N = 5). The charge density is rather insensitive to the 
SO coupling, a. However, we show that the SO interac- 
tion couples the spin polarization to the orbital motion 
determining the spin properties of the GS and the first 
excited states in a surprising way. Indeed , by increasing 
a, the expectation value of the spin density of the GS, 
which was originally oriented in the z— direction, acquires 
a component in the dot plane, because the minority spin 
density is increased and pushed from the center of the 
dot outward. Moreover the combined effect of U and a 
deforms substantially the spin density of the SKD. A 
sharp minority spin polarization is present close to the 
dot center. The reversal of the spin polarization at the 
origin in the SKD state w.r. to the FSP GS leads to an 
extra node in the spin density. 

This situation is reminiscent of the case of the Quan- 
tum Hall Ferromagnet (QHF)2& close to filling one. In 
that case, a true magnetic ordering is achieved, which 
is characterized by full spin polarization in the GS and 
by a topologically constrained first excited state, the 
Skyrmion (SK) state, with reversal of the spin at the ori- 
gin, first studied in the 0(3) non linear a model (NLaM) 
in 2— d dimensionsS&Si. In Section V, we elaborate on the 
analogies and differences between the FSP dot and the 
QHF. In the QHF a topological quantization of charge 
occurs. By contrast the SKD state has no topological 
features, because the geometrical compactification pro- 
cedure described in Subsection V.A cannot take place. 
We refer to the SKD state as a "spin exciton" because 
there is some piling up of the charge at the center of the 
dot w.r. to the GS, together with the reversal of the spin 
polarization there. 

The paper is organized as follows: 

In Sect. II we report our results of numerical diagonal- 
ization with the Lanczos algorithm for a dot with five 
electrons (N = 5) close to the FSP state, in the absence 
of SO coupling. 

In Sect. Ill we derive the SO matrix elements in the 
2 — d harmonic oscillator basis and discuss how the SO 
coupling modifies the lowest lying energy levels. 

In Sect. IV we show the spin and charge density of the 



lowest lying J z multiplet when the FSP state is achieved. 

In Sect. V we summarize the features of the SK state in 
the QHF and compare these with the ones of the SKD 
state in a dot with SO coupling. 

A brief summary and some conclusions are outlined in 
Sect. VI. 

There is evidence of skyrmion excitations in Ga As 2 — d 
electron gas systems close to filling one by magnetoab- 
sorption spectroscopy^. A sharp absorption line could 
be found in exciting dots to the SKD state, by trans- 
ferral of energy and angular momentum with circularly 
polarized light. This amounts to adding a spin exciton 
to the dot. 



II. FSP STATE AND DOT RECONSTRUCTION 

We consider N = 5 electrons confined in two dimen- 
sions (spanned by the (p, tp) coordinates) by a parabolic 
potential of characteristic frequency u>d- This is a model 
for an isolated disk shaped QD. A magnetic field B or- 
thogonal to the disk is measured in units of huj c (meV) , 
where oj c is the cyclotron frequency. In the absence of 
spin orbit coupling, the single particle states <f) nm are the 
eigenfunctions of the 2 — d harmonic oscillator with fre- 



quency lo — <J uj 2 , + -j-. They are labeled by n, m (with 
n G (0, 1, 2, 3, ...) and m £ (— n, — n + 2, n — 2, n)). m 
is the angular momentum in the z direction: 
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Here L" (t) (with t = p 2 /l 2 ) is the generalized Laguerre 
polynomial with n > I = yjh/m*ui is the char- 
acteristic length due to the the lateral geometrical con- 
finement in the dot inclusive of the B field effects and 
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is a normalization factor. 



The corresponding single particle energy levels are 

jn 

e n ,m = {n+ l)huj - —fiuj c . 



(2) 



In the absence of both interaction and magnetic field, the 
lowest lying single particle states are occupied with the 
minimum spin. The GS Slater determinant is sketched 
pictorially in Fig^,, where energy is intended on the ver- 
tical axis. Each box represents a single particle state la- 
beled by n, m and arrows represent electron occupancy 
with spin projection along the quantization axis. Wc 
have performed exact diagonalization of this system in- 
cluding Coulomb interaction between the electrons. The 
matrix elements of the unscreened Coulomb interaction 
use the single particle basis set, up to 28 orbitals. They 
can be calculated analytically and are parametrized by 
the strength of the interaction U . Our calculation is lim- 
ited to very small particle numbers [N < 7), because the 
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a) GS at B = 



b) GS at B = B' 



FIG. 1: Slater determinants quoted in the text are depicted. 
Quantum numbers are N = 5, S — 1/2 for the state at B = 
[a)] and S — 5/2 for the state at B = B* , the magnetic field 
value at which the maximum of S is achieved [6)1. 



truncation of the Hilbert space influences the results for 
larger N. However, our convergency checks show that 
the numerical errors proliferate only at higher energies. 
In particular they affect the reliability of the level spin 
degeneracy. In any case, numerical errors are quite small 
if is N = 5. In Fig|2] [left panels], we show the lowest 
lying total energy levels at fixed angular momentum M, 
versus M, for U = 13meV. Magnetic field is B = 5meV 
[top], B = B* = 7meV [middle], B = U.SmeV [bottom]. 
At each M, the spin degeneracy is marked by dashes of 
different length: short dashes for S — 1/2 (doubly de- 
generate level), medium dashes for S = 3/2 (fourfold 
degeneracy) and long dashes for S — 5/2 (sixfold degen- 
eracy). On the r.h.s. of the picture the radial charge 
density of the corresponding GS is plotted vs distance r 
from the dot center. FigEl ([left panels]) shows the cross- 
ing of levels with increasing B. Electron-electron corre- 
lations imply that when M increases, S also increases. 
At B = B* = 7 meV the spin S reaches its maximum 
value S — N/2. The largest contribution to the GS wave- 
function is given by the Slater determinant depicted in 
FigOJ)) corresponding to M — ^ _1 to — 10. We con- 
centrate on the state at B = B*,the FSP GS. This cor- 
responds to the "maximum density droplet" state dis- 
cussed in the literature^*. Qualitatively we can say that 
at B = B* the dot attains its smallest radius. As can be 
seen from the GS charge density, further increase of B 
leads to the so called reconstruction of the charge den- 
sity of the dot. For B > B* , the M of the GS increases 
further, but S is no longer at its maximum. In the bot- 
tom panel of FigE] it is shown that at B = 11.5 meV 
the GS energy is now achieved for M = 13 with a dou- 
blet (S = 1/2) state. The corresponding charge density 
of the dot, as depicted on the r.h.s, is strongly modified 
close to the edge 15 : it displays a node followed by an 
extra non zero annulus at a larger distance. In view of 
the fact that our expansion of the wavefunction only in- 
cludes rotationally invariant components, the breaking of 
the azimuthal symmetry is impossible. By contrast this 
is found to occur in density functional calculations and 
the corresponding GS is referred to as the de Chamon- 
Wen phasei*. The GS at B = B* can be compared 
with a FSP quantum Hall state of an extended disk in 
the absence of lateral confinement (Quantum Hall Ferro- 
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FIG. 2: Energy levels without SO coupling for the dot with 
N — 5 electrons at U = 13 meV and u)d = meV. Magnetic 
field values are: (in units of hui c ): B — 5meV [top],£> = B* — 
7 meV [middle], B = 11.5 meV [bottom]. The total M is on 
the x axis. The levels are drawn with short, medium or long 
dashes, depending on the total spin : S = 1/2, 3/2, 5/2 



magnet (QHF) at filling one). Fig. QJd) recalls the occu- 
pancy of the lowest Landau level (LLL) up to a maximum 
to = N — 1, except for the fact that in our case the single 
particle levels corresponding to the LLL are not all de- 
generate in energy. In the language of the quantum Hall 
effect the unperturbed levels are: 



= [2v- 



TO 

l)huj - —hue 



(3) 



where v — (n — |to|)/2. LLL is for v = and to > 0. 
The Slater determinant of Fig^) has a charge density 
which is flat as a function of the radius r, up to the disk 
edge, where it rapidly falls down to zero. In our case this 
feature is lost because of the presence of U, together with 
the fact that the number of electrons is small. We will 
better discuss the comparison of the FSP GS with the 
QHF in Section V. 



III. INCLUSION OF SPIN-ORBIT 

We now add the spin orbit interaction to the FSP GS 
at B = B* . This can be tuned by applying an electric 
field E in the z direction, which couples to the spin of 
the electrons in the dot with a term 21 : 



(z x p) 



a . 



(4) 



Here a are the Pauli matrices, a is the spin-orbit coupling 
parameter which is proportional to the electric field, a 
will be measured in units of meV ■ A. We now rewrite the 
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spin-orbit coupling term in a second quantized form. We by c nma , cjj mtr and we get: 
denote the fermion operators associated to <f) nm of Eq. Q 
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< n m 
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(5) 



The integration over the azimuthal angle ip can be done 
analytically. This shows that the Hamiltonian can be 
rewritten in the following way: 

nn' m 

An 1 m-l,™^' m— l^ nm i-^ (^) 



with A 



n'm'nm 



f°° dm 
S m '+i,m J dtR n ,\ m ,\(i)(2\/i— + -j=)R n \ m \(t) 



f°° d'm' 
$m'-i,m J dtR n ,\ m >\(t)(2— Vt + —jz)R n \ m \ (t) 

Here -B„ m ,„' m _i = A* n , m _ l nm , what implies that the 



hamiltonian is hermitian. It is clear that while s z and m 
are no longer separately conserved, their sum j z = s z +m 
( with j z half integer ) is a good quantum number. We 
will denote the single particle basis that diagonalizes the 
SO term by with (3 — p,m. The label (3 takes two 
possible values, say p, q and allows for conservation of the 
number of degrees of freedom. 

The SO interaction lifts the spin degeneracy. In FigEl 
we show the splitting of the multiplet with N = 5, 
S = 5/2, M = 10 at B = B* = 7 meV and U = 13 meV 
vs the strength of the SO coupling a. The strength of U 
is responsable not only for the fact that the GS belongs 
to this multiplet, but also for the ordering in energy of 
the sequence: J z = 25/2,23/2,21/2,19/2,17/2,15/2 ( 
from bottom to top). At small U values the sequence 
is J z = 15/2,17/2,19/2,21/2,25/2,23/2, as shown in 
Fig^] With increasing of U, some level crossings oc- 
cur. The ordering at three different values of U is mag- 
nified in the bottom panels of FigQ] The case with 
U = 13 meV is shown in the bottom right panel of 
Fig. 0] The lowest state in energy is for J z = 25/2 fol- 
lowed by J z — 23/2, 21/2 (almost degenerate with 15/2), 
15/2, 19/2, 17/2. At U = 13meV a sizeable gap is formed 
between the J z — 25/2 GS and the first excited state 
J z = 23/2 (SKD). The other states of the multiplet are 
bunched together at higher energy. In this Section we 
focus on the U — 13 meV case and discuss the charge 
density and the spin polarization density of the GS at 
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FIG. 3: Splitting of the lowest lying multiplet for iV = 5,5 = 
5/2 and M — 10 vs strength a of the SO interaction, at 
B = B* = 7 meV, U = 13 meV and uj d = 5 meV. The levels 
are labeled by J z . 



J z = 25/2. The other states of the multiplet will be 
analyzed in Section IV. 

As it appears from Fig. [S] [top panel], the charge den- 
sity of the GS is only mildly changed when we increase 
the SO coupling. By contrast, the spin density is quite 
sensitive to the addition of SO, up to saturation. Now 
the z— component of the total spin is no longer a good 
quantum number and some admixture with down spin 
electrons appears. Indeed the role of the Rashba term is 
to rotate the average electron spin. In particular, down 
spin electrons are pushed away from the center of the 
dot, giving rise to the spin density components S z (f) (or- 
thogonal to the dot plane), and S r {f) (in the plane of the 
dot), which are plotted in Fig[S] [middle and bottom pan- 
els, respectively]. It is remarkable that the spin density 
S z (r) changes sign at the edge of the dot for large SO 
coupling. This is confirmed by a plot of the occupation 
numbers n nma = (G5'|cJ lmo .c nm(J |G5') with n = m. They 
are shown in Figl^Jfor both N = 4 and N = 5 for com- 
parison. Of course, the change of N would also imply an 
effective change of the confinement potential ujd (what 
we do not do). However, all what we want to show here 
is that our findings depend on the strength of B only, 
and not on the number of electrons being even or odd. 
A similar feature occurs in the de Chamon-Wen phase, 
in the absence of SO: when crossing the edges, the spins 
tilt away from their bulk directional. 

The reversal of the spins in the tail at the dot bound- 
ary is a peculiarity of the Rashba interaction, but the 
spin/charge density is very small there and does not in- 
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FIG. 4: Energy levels with B — 7 meV, u)d — 5 meV, and 
a = 100 meV ■ A, for different (7 values. In the upper panel 
we show the crossings that allow the FSP polarized state to 
be the ground state when U is large. Ordering of the levels is 
magnified in the bottom panels for three different U values. 
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FIG. 6: Occupation numbers n n=mim , CT in the GS with N = 
4(5) electrons (left(right)) without SO (top) and with SO 
(a — 100 meV ■ A) (bottom). Other parameter values are 
B — 7 meV, U — 13 meV, uid = 5 meV . White bars refers to 
spin up, grey bars refer to spin down. The FSP GS of the dot 
with N — 4(5) electrons has total spin S = 2(5/2) and the 
2— component of the total angular momentum J z — 8(25/2). 



fluence the dot properties. 



IV. SPIN AND CHARGE DENSITY IN THE 
MULTIPLET 5 = 5/2, M = 10 
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FIG. 5: Charge density, azimuthal spin density S z , in plane 
spin density S r , in the radial direction, in the GS (J = 25/2) 
at various SO couplings: a = 5, 100, 250 meV ■ A. Here B = 
7 meV , U = 13 meV and uJd = 5meV 



In the previous Section we have shown that at B = B* 
the GS with N = 5 electrons belongs to the S = 5/2, 
M = 10 multiplet. The SO coupling lifts its degeneracy 
as shown in FigEl The size of U strongly influences the 
energy of each state, by producing crossings of levels. 
At U = 13 meV the lowest lying states with increasing 
energies are (see FigQJright bottom panel]): 



\GS>= \N = 5; 5 = 5/2, J* 
spin polarized GS . 



25/2 >: this is the fully 



\SKD >= \N 
exciton'. 



5; S = 3/2, J z = 23/2 >: the 'spin 



\b >= \N = 5; S = 3/2, J z 
higher in energy w.r. to 



= 15/2 >: this is a state 
SKD >. 



This ordering of energy levels is again a consequence of 
Hund's rule: lowest energy is for J z = L z + \S Z \, higher 
energy is for J z — L Z ~\S Z \. Besides affecting the energy 
of the states, the effect of U is to enhance the trans- 
fer of weight from the majority ("up") to the minority 
("down") spin population. This is shown in Fig0 where 
the occupation numbers n n=m ^ m . a are reported for the 
states \GS >, \SKD > and |6 >'for U = [left panels] 
and U = 13 meV [right panels], respectively. A striking 
feature characterizes the spin densities of these states (see 
Figd [SJ: the dominant spin density is reversed in the 
\b > state, w.r. to the \GS >. The state \SKD >, which 
is the first excited state, interpolates between the two. 
Spin occupancy is not significantly modified for larger r. 
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FIG. 7: The occupation numbers n n = m , mjCT in the state at 
J z = 15/2,23/2,25/2 for small U [left], and large U [right]. 
White bars refers to spin up, grey bars refers to spin down. We 
stress that at U — the ordering of the levels, corresponding 
to the three panels on the left is changed w.r. to the ones on 
the right (U = 13meV). (see Fig. [I][bottom panels]). 




While at U = the flipping of the spin at the origin w.r. 
to the GS is full, in the interacting case some up-spin 
is left at the center. This allows for a smoother radial 
dependence of the spin and charge density expectation 
values. Eventually, this is the reason why this state turns 
out to be the lowest excited state in the FSP system. 

In FigEl w e show the charge and spin densities of the 
complete multiplct at a = lQOmeV ■ A, U = ISmeV and 
B = B* . The situation is quite peculiar: by looking at 
< S z > [middle panel], we see that the GS has an up 
spin density everywhere in the dot, except for a tiny lit- 
tle reversed tail at the boundary. By contrast, the state 
J z = 15/2 has a down spin density at any r. Intermedi- 
ate between the two, the SKD state displays a reversed 
spin at the center of the dot but the spin polarization 
changes into up when approaching the edge, to restore 
the spin density of the 25/2 state. There is a node in 
the middle. The other states (17/2, 19/2 and 21/2) are 
rather featureless and they do not share these features. 
The trend is confirmed by looking at the projection of 
the spin density in the plane of the dot S r — f ■ S (see 
FiglSJbottom panel]). This is the complementary infor- 
mation w.r. to S z (r). When S r (r) in strongly non zero, 
then S z (r) is heavily reduced. 

An analogous interpolation occurs for the charge den- 
sity. There is a piling up of the charge at the origin (see 
FiglSftop panel]), corresponding to a locally dominant 
down spin density. The [SKD > state is a collective ex- 
citation of the QD, which we call a "spin exciton". In 
the next Section we show that the spin exciton recalls 
the first excited state of a QHF with some important 
differences, though. 



FIG. 8: (color on line) Charge density, azimuthal spin 
density S z , in plane spin density S r , in the radial di- 
rection, at various J z . From bottom to top: J z = 
25/2(GS), 23/2(SKD), 21/2, 19/2, 17/2, 15/2. Other param- 
eters are a — 100 meV ■ A, U = 13 meV and ujd = 5 meV 



V. COMPARISON BETWEEN THE DOT AND 
A QHF DISK 



The case of the dot in the FSP state can be compared 
with that of a disk shaped Quantum Hall Fcrromagnet 
at filling one. The comparison is in order, because the 
physics of the dot turns into that of a quantum Hall 
disk by increasing the magnetic field, as long as the ratio 
LUd/iUc — * 0. Of course, while the infinite quantum Hall 
system is marked by a phase transition to the spin polar- 
ized state, the dot, being a system with a finite number of 
particles, undergoes a crossover to the FSP state which 
is not a broken symmetry state. This is confirmed by the 
presence of the tiny minority spin tail at the edge of the 
dot. 

In the Subsection V.A, we recall some properties of 
the Hartree Fock description of the GS and first excited 
state of the QHF, which applies to filling close to (but 
less than) one. 

Similarly,somc analytical approximations leading to a 
simplified H — _F-like approach for the dot with SO cou- 
pling will be discussed in Subsection V.B to highlight the 
analogies between the two systems. 
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A. Quantum Hall Ferromagnet 

In describing the QH state on a disk it is customary 
to label one particle states with v = (n — \m\)/2 and 
m, a, corresponding to the eigenvalues e„ m CT given in 
Eq.@. The LLL includes the wavefunctions 4> nm given 
by Eq.QJ for v — and m > 0. In this case all Laguerre 

polynomials Lq (t) = 1. If there is no confinement po- 
tential (bJd = 0), all eom are degenerate. We rename the 
LLL wavefunctions fomXo- (here x<? denotes the spin 1/2 
wavefunction) and associate the single particle fermion 
operators a u =o ma to them. In the QHF at filling one, 
the LLL subband with, say, spin up, is fully occupied: 
the GS is a fully polarized spin state: 



QHF.O 



n 



*0mt 



0<m<N-l 



(7) 



Here \vac) is the vacuum state. The lowest lying branch 
of excitations of the QHF are spin waves. These involve 
electrons in the down spin LLL subband and holes in the 
up spin LLL subband. 

It was pointed out long ago^ that, if the filling is 
slightly less than one, the first excited state can be a 
very special collective excitation with S < N/2 and an 
extra node in the spin density. The spin polarization is 
reversed at the center, but gradually heals to the domi- 
nant spin background over a distance of many magnetic 
lengths (SK state). This excitation can be traced back 
to the skyrmion , the topological excitation of the 0(3) 
NLaM in 2 — dM. A disk of infinite radius in coordinate 
space can be compactified to a sphere S 2 in 1Z 3 having 
the origin in the south pole and the point at infinity in the 
north pole. A similar compactification can be performed 
in the order parameter configurational space. An uniform 
magnetization "up" is represented by a vector pointing 
to the north pole everywhere on S 2 . The skyrmion is 
a finite action configuration on S 2 , satisfying the classi- 
cal eq.s of motion for the magnetization of the NLaM, 
conserving J = S + M and belonging to a non trivial 
homotopy class. If the topological charge is Q = 1, the 
shape of the magnetization field is s(r) = f, where f is 
the normal to S 2 at each point. Q is the flux of s(r) 
through the sphere of unit radius. The spin polarization 
is "down" at the south pole and turns over continuously 
in space, until it reaches "up" at the north pole. That is, 
the spin polarization is flipped at the origin of the disk 
w.r.to the GS and turns smoothly over away from it in 
the radial direction. 

Within Hartree-FockS, the Slater determinant \S,K) 
that describes this state conserves total J z . To con- 
struct it, a canonical transformation is performed on the 
fermion operators: 

.1 



P-t = a QOl 



j e (-, ...oo) 



-Vj a 0i _i T + uj a 0j+ i i , j € (-, ...oo) ,(8) 



Normalization requires that |it,- 



1. Note that 



the operator p_ i still belongs to the LLL as it destroys a 

particle in the /i/ = o, m =oXl state. We denote by /J' 9 the 
single particle orbitals corresponding to the operators of 
Eq.© and we use them in Appendix A. 

The generic Slater determinant built by means of these 
operators is: 



OO p q 

\S,K)=X\(pU) n ^^\ V ac) 



(9) 



are the occupation numbers of the single particle 



states (n 



p 



--< p)pj >,n g j 



< >)j with E 



if> n 3 







N . The state of Eq.® is labeled by the total spin S and 
by K . S z is no longer a good quantum number and is 
substituted by 



K = S 



N/2 

-Y 



(n? 



(10) 



The state of Eq.© with S = N/2,K = is the FSP 
QHF ground state of Eq.©, if the only non zero oc- 
cupation numbers are n? = 1 for j G (i, 7V/2) with 
uj = 1 for j E (^, N/2). This state corresponds to the 
FPS GS of Figt!J>) for the QD case. 
For the hard core model the HF equations can be 
solved analytically 2 ^. The lowest lying skyrmion state 
is l-W/2,1), with 



M 2 = i-H 2 



e 



e + c? + |) 



(ii) 



leading to the spin density s(r) defined in terms of the 
arbitrary length scale £ ( r 2 = x 2 + y 2 )24 (see Appendix): 



Sxif) 



2x£, 



2 +e 



s z (r) = 



e 



(12) 

The ± refer to the sign of the topological charge Q = ±1. 
In the real QHF the length £ is governed by the relative 
strength of the Zeeman and the Coulomb energies. 



B. Dot with spin-orbit coupling 

In this Subsection we give arguments supporting our 
claim that the state SKD of Section IV corresponds to 
the state \N/2, 1) in the QHF limit, that is in the limit 
of zero confinement potential and filling one. Indeed, the 
radial distribution of the spin density of the state SKD 
recalls the one of Eq. (|12fl except for a very shallow tail at 
the boundary. Away from the center the spin polarization 
of the SKD state lines up gradually with the one of the 
GS as it happens for the case of the Skyrmion. As in 
the SK case, S z (r) has an extra node at r = f. In the 
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presence of SO, the length scale £ is no longer arbitrary, 
but is fixed by the strength of the SO coupling. 

In the case of the QHF on a disk, both rotations in real 
space around the z— axis and rotations in spin space are 
good symmetries, so that M as well as S z are conserved. 
This implies that an allowed SK— like excited state of the 
real system has to be obtained by projecting the state of 
Eq.@ onto the subspace of definite M and S z . This is 
not necessary in the QD with SO interaction, because 
the SO hamiltonian term only conserves J z as the state 
\S, K) does. In the following we show that a simplified 
H — F'-like approach for the dot case with SO coupling 
shows features similar to the ones described by Eq.s JSJ), 
jnj and Eq.s ©, JEJ. Let us first discuss SO coupling 
in the dot at U = 0. The vector space required to diago- 
nalize the SO coupling and to obtain the eigenfunctions 
Wj exceeds the LLL space enormously (in practice we 
always use the basis of Eq.QJ and never calculate the 
Wj^ 's explicitly) . As a simple analytical approximation, 
we can restrict ourselves to the LLL for sake of simplicity. 
We have checked numerically that this approximation is 
largely satisfactory away from the level crossings. In this 
case, diagonalization of the SO interaction factorizes the 
problem into a collection of 2 x 2 matrices. What the 
SO does is to mix single particle states with different m 
and opposite spins in the way that the transformation 
of Eq.JHJ) shows. Indeed, j z (j z = j in the following) is 
conserved. Within the LLL, two (to, a) values contribute 
to each half integer j: (to, |) and (to + 1, J,). The unper- 
turbed energy levels involved, eo m and eo m +i, are given 
by Eq.©. Let the offdiagonal matrix element including 
the SO coupling be a. Then the eigenvalues are: 



= -(e 0m + e 0m+ i)±y T +« 2 , (13) 

where S = eo m +i — eom = — ^ c /2. The di- 
agonalization implies a rotation in the 2-vector space 
{/omXT> /om+iX|}of an angle 7 given by tan 27 = -2a/ 5. 
The single particle states obtained in this way coincide 
with ff q defined after Eq.©. The mixing of the two 
states (to, I) and (to + 1,1) is j— independent, within 
our approximations, because 5 is. This implies that the 
rotation angle 7 keeps roughly constant in the radial di- 
rection, because average radial distribution of an elec- 
tron of angular momentum j is ~ ly/j + 1. We can now 
construct the Slater determinants representing the states 
lower lying in energy. The states corresponding to the 
ones obtained numerically in Section IV are depicted in 
Fig.@. In analogy to Fig.Q we use boxes to allocate 
electrons. Each box is cut into a lower and an upper tri- 
angle w.r.to the diagonal, corresponding to the q and the 
p state of a given j z , respectively. A dot marks which 
of the orbitals is occupied. We have analyzed the Slater 
determinants which contribute mostly to the states ob- 
tained at the end of the Lanczos procedure, giving the 
average occupation numbers of Fig.Q. Their largest 
components indeed contain the determinants shown in 




- 1/2 - 1/2 - 1/2 - 1/2 



GS SKD a) b) 

FIG. 9: Slater determinants quoted in the text with the same 
labels. Quantum numbers are N = 5, S = 5/2 and J z . 
Upper/ lower triangle refer to single particle states labeled 
by j z and p/q. Other possible quantum numbers do not ap- 
pear. The dots mark occupied states. Configuration labeled 
as a) is involved in a state belonging to a much higher energy. 

Fig.©. 

There is close similarity with the results of Subsection 
V.A. However a relevant difference can be immediately 
recognized. While the skyrmion shows a very smooth 
tilting of the spin orientation with increasing distance 
from the center of the disk (see Ea. (|ll|l % ). the rotation 
angle 7 for the dot is uniform in the radial direction. 
This feature is partly compensated by the addition of the 
Coulomb repulsion. Indeed, U 7^ predominantly affects 
the occupations close to the center of the dot disk, while 
its influence fades out away at larger distances. This fact 
introduces a radial variation of the tilting of the spin 
polarization. According to Ea. H12l) . the skyrmion has a 
linear variation with radial distance of S r , close to the 
origin. By contrast, our numerical results reported in 
Fig. (0 [bottom panel ]) show a quadratic increase at small 
r's. 

The role of U is quite substantial, by locating the en- 
ergy of the SKD state intermediate between those of the 
GS and of the b state. Needless to say, another relevant 
and obvious difference between SKD and SK is the ab- 
sence of any conserved topological charge in the dot. In 
the QHF the conservation of Q is implemented on sym- 
metry grounds, by mapping the QH disk onto a sphere. 
This mapping cannot be extended to the dot, because, as 
seen from Fig.JSJ), the direction of the magnetization at 
the boundary is not unique. Magnetization is not defined 
at r — ► 00: the point at infinity is a singular point in the 
magnetization configurational space. 



VI. CONCLUSIONS 

In a disk shaped quantum dot with few electrons, in- 
teractions drive the system to a fully spin polarized state 
with S = N/2, in the presence of a magnetic field B = B* 
orthogonal to the dot. At B > B* , the total spin is again 
drastically reduced and the charge density reconstructs 
at the disk boundary. We have reported on exact diag- 
onalization results of a QD with N = 5 electrons and 
studied the effect of SO coupling possibly due to an ex- 
ternal electric field orthogonal to the dot disk. There 
are analogies between the dot state at B = B* and the 
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Quantum Hall Ferromagnet at filling one. We require a 
sizeable interaction strength U to stabilize the FSP GS. 
When the SO coupling is increased, level crossings occur 
in the splitted 5 = 5/2 multiplct, until the state with 
maximum J z = M + S z = 25/2 becomes the GS. The 
first excited state (SKD state) has J z = 23/2. When 
compared to the GS, the SKD state has some charge 
tranferred to the dot center and a very peculiar spin tex- 
ture. Indeed, the z— component of the spin density at 
the center of the dot is opposite to the one of the GS 
and rotates continuously over away from the center, by 
acquiring the same profile as the one of the GS at the 
dot boundary. This winding requires an extra node in 
the spin density, which is absent in the other multiplct 
states. According to these properties, the SKD state 
can be viewed as carrying one spin exciton. Both our 
numerical results of Section III-IV, and our approximate 
analytical speculation of Section V show how essential 
the combined role of the SO coupling and of the e — e 
interaction is in stabilizing this state. Our calculation 
parametrizes the interaction strength U, but it ignores 
the screening of the e — e interaction altogether. This 
should be reconsidered in view of the fact that vertical 
QD's are separated on the top and the bottom from the 
contact metals by barriers with a typical width of 70A. 
Even for N = 5 this is smaller than the inter-electron 
spacing 30 . However, we believe that the exciton state is 
robust when the screening is included. Indeed, the flip- 
ping of the spin is concentrated at the center of the dot 
and is governed by the e — e interaction at short range, 
which should be largely insensitive of screening effects. 
The SKD state recalls the skyrmion excitation which 
takes place in a disk shaped QHF at filling one. The 
statement could be puzzling, in view of the fact that the 
SO coupling is essential to the SKD state, but it is never 
invoked when discussing Quantum Hall properties. How- 
ever field theory models (NLaM) use the conservation 
of J to prove the existence of the skyrmion state. In a 
real isolated QH disk M, S z would keep finite and sepa- 
rately conserved. In this case only the component of the 
SK state that conserves given values of M, S z would be 
present in the excitation spectrum. Nonetheless the dif- 
ference is washed out in the limit of an infinite disk size. 
This is the continuous limit which leads to the NLaM. 
In the case of the dot, the compactification of both the 
coordinate space and the magnetization space cannot be 
performed because the direction of the magnetization is 
not defined at r — > oo. Therefore no state can be con- 
structed that conserves J z only, without conserving M 
and S z separately. The spin orbit coupling opens up this 
possibility. However, no topological charge can be asso- 
ciated to the SKD state in the dot. 
Our calculation shows that for realistic values of the dot 
confining potential (w<j = bmeV), of the Coulomb inter- 
action strength (U = 13meV) and of the SO coupling 
a ~ lOOmeVML, the FSP GS and the SKD state arc 
well spaced levels. The other levels of the multiplet ap- 
pear at higher energies and are rather close to each other. 



This means that, at B = B* , the dot opens a sizeable 
spin gap between the GS and the SKD state, that can 
be tuned with an applied gate. This spin gap cannot be 
washed out by thermal fluctuations, if the temperature 
is low enough (~ 50mK). The gap can be probed by op- 
tically pumped NMR as in quantum wells 22 . Spin-lattice 
relaxation of 71 Ga nuclear spins in the dot, driven by the 
hyperfine coupling to the dot electrons should be very 
much reduced, thus leading to a large T\ . 
The extremely low spin relaxation expected for this ex- 
citation, could allow for a coherent manipulation of the 
spin exciton using terahertz radiation 3 ^ 3 ^ 7 .. In general, 
we believe that the system studied here can be relevant to 
the coherent manipulation of QD states. This is appeal- 
ing in view of quantum information processing 32,33 34 . 
Indeed, a spectrum like the one calculated in this work 
should produce sharp optical absorption lines. Photolu- 
minescence induced by a pump and probe laser technique 
has been studied in disk shaped In Ga As QD's with ev- 
idence for Rabi oscillations 3 ^. In our case, because of the 
presence of B, a circularly polarized pulse of one single 
chirality can excite the spin exciton discussed here. 
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APPENDIX A: QHF SPIN DENSITY 

In this appendix we show that the state \N/2, 1) given 
by Eq.@ with Uj given by Eq.JTT} leads to the skyrmion 
spin density of Eq. (|I2|I . 

The wavefunctions for the QH disk associated to the 
operators a,Q ma are given in Eq.Q. In the LLL [y — 
(n — \m\)/2 = 0) all Laguerre polynomials L™ (t) = 1. 
To construct the field operator, we associate a spino- 
rial wavefunction f^ q {r) to the operator pj/qj following 
Eq.©: 

W ,6( i-~> (Al) 

and analogously for /?. We take Uj and Vj real. The field 
operator is: 

oo 

4(f) = E (//-!(*%•-! + /'(*%■) • (A2) 

l = i 

The spin density operator is s*(r) = 3?e |^(r)cn/>(r) |, 

to be evaluated on the state \N/2, 1). Let us consider 
s x (r) first. The term including the jij operators does not 
contribute, because all the f p orbitals are unoccupied in 
the state \N/2, 1), except for j = — i. On the other hand 
this term, does not appear, because U_i -v_i =1-0 = 0. 
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The contribution to s x (r) given by the qj operators is: Because the maximum of the first factor occurs for j+| ~ 

r 2 we evaluate the denominator of the second factor by 

E( w i/oj_i(r) u i/o-+i(r>) ( 1 \ ( u ifoj \ substituting j + | — > r 2 , what allows us to perform the 
~ 3 2 \ 1 y I v jfoj+^(^) / sum ex P nc itly By noting that rcos</? = a; we obtain 

J = 2 " s x (f) as given by Ea. (|12jl . A similar calculation applies 

00 for s y (r). In the case of s z (r), the extra factor is v 2 — it 2 = 

= zl 2u i v jfaj-±(?)foj+± ( f ) ■ ( A 3) [(j + i) _ + i) + and ^ disappears. Using the 

i=3 same approximations as above, we obtain the result of 

_ Eq.(0. 
Using Eq.JTTJ) we get (r = (r, <£>)): 

r i-| r j+i (j + i)3 _ r2 | 




(j-l)! € 2 + a + ' 



(A4) 
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